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A minimal model is proposed for the perovskite manganese oxides showing the strongly in- 
coherent charge dynamics with a suppressed Drude weight in the ferromagnetic and metallic 
phase near the insulator. We investigate a generalized double-exchange model including three 
elements; the orbital degeneracy of Eg conduction bands, the Coulomb interaction and fluctuat- 
ing Jahn- Teller distortions. We demonstrate that Lanczos diagonalization calculations combined 
with Monte Carlo sampling of the largely fluctuating lattice distortions result in the optical con- 
ductivity which quantitatively accounts for the experimental indications. It is found that all the 
three elements are indispensable to understand the charge dynamics in these compounds. 

KEYWORDS: perovskite manganites, double-exchange model, Coulomb interaction, orbital degeneracy, Jahn- 
Teller distortion, metal-insulator transition, ferromagnetic metal, optical conductivity, Drude weight, 
incoherent charge dynamics, exact diagonalization, strongly correlated electron system 
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As an open problem of perovskite manganese oxides 
Ri-^cA^cMnOs (i?=rare earth, A=Ca, Sr, Ba WhRUm 
beyond the double-exchange (DE) mechanism;BBa'B'EP 
strongly incoherent charge dynamics in the optical con- 
ductivity (j{uj) of a typical compound, Lai^^Sr^MnOs, 
in the ferromagnetic and metallic (EMl phase near the 
insulator (x=0.175, 0.18) is knownE'ElEF The coherent 
component of cr(a;), namely the Drude weight, at low 
temperatures is suppressed; its ratio to the effective car- 
rier density is about 0.2 and the nominal mass enhance- 
ment is very large (m*/m=50-80). The two-dimensional 
(2D) system of Lai.2Sri_JVIn207 also reveals the inco- 
herent charge dynamics Jl3) Such an incoherent feature 
cannot originate from the spin degrees of freedom under 
the spin polarization in the FM phase. 

Possible elements as candidates contributing to the 
incoherence are orbital degeneracy, electron correlation 
and dynamical Jahn- Teller (JT) distortion, which are 
neglected in the simple DE model. A question we are 
faced with is what is the minimal model for the mangan- 
ites showing such a strong incoherence. Previous stud- 
ies suggest that a partial account of the three elements 
cannot reproduce the dominance of incoherence quanti- 
tatively, (i) Only the eg-orbital degeneracyliiP leads to 
an incoherence comparable to the coherent part; the ra- 
tio of the Drude weight to the total weight is ^0.46. (ii) 
The e--£ixbital degeneracy combined with the JT distor- 
tionalj^H gives the ratio ~0.7. (iii) For a model contain- 
ing both the eg-orbital degeneracy and the electron corre- 
lation, a broad incoherent structure of cr{uj) in the energy 
region below the charge gap for the slightly hole-doppjl 
case is obtained. However, the mean-field treatmentEfP 
gives the absent Drude weight. Numerical calculations in 
2D at substantially high temperature conclude the ratio 
of the Drude weight to the effective carrier density to be 
~0.4, in spite that the coherence is quite suppressed due 
to the thermal fluctuations. Ej^ Here the effective car- 



rier density is defined as the Drude weight added by the 
incoherent weight within the Wjpi' Hubbard band. Nu- 
merical calculations at T = 0ji3'll3^ where the quantum 
fluctuations are fully taken into account, give the ratio 
~0.45 under the strong Coulomb interaction, even when 
the chemical potential difference between two orbitals in 
a 2D plane is controlled to let the orbital polarization 
(OP) vanish so that the three-dimensional (3D) system 
is simulated. Thus, if only a part of the three elements 
is taken into account, the Drude and incoherent weights 
are at most comparable and cannot reproduce the domi- 
nance of incoherence observed in charge dynamics of the 
experiments. 

In such circumstances, a model containing all the three 
elements, the e^-orbital degeneracy, the electron correla- 
tion and the JT distortion with the staggered and long- 
ranged pattern of the displacement of the oxygens sur- 
rounding the Mn atoms has been investigated. A coop- 
erative effect of these elements are known to play an es- 
sential role for tJaeJundoped compounds showing an insu- 
lating behaviorlla'ta* The results for the doped metalEl^ 
suggest that this staggered distortion under the strong 
correlation makes the charge dynamics more incoherent 
because the distortion enhances the orbital fluctuation 
obstructing the motion of electrons. In the experiments, 
on the contrary, no long-ranged orders of JT distortion 
are observed at the doping where the strong incoherence 
appears, although short-ranged JT correlations are ex- 
pected to persist. To reach full understanding of the ori- 
gin of the incoherence, it is required to examine whether 
the charge dynamics becomes as incoherent as observed 
in the experiments when we treat such short-ranged JT 
correlations together with the orbital degeneracy and the 
electron correlation. 

When the displacement of the oxygens is considered, 
one should be careful about a zero-point energy (ZPE) 
which comes from the quantum nature of the lattice. 
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One can decompose the displacement of the oxygens 
into the normal modes of the lattice, among which we 
consider the Q2 and Qs modes coupling with the elec- 
trons in the Cg orbitals through the JT coupling. Conse- 
quently, one has the phonon part in the Hamiltonian to 
be T.^[-J^{^ + ^) + k{(l^^ + ql^)] whcrc mo is the 
mass of an oxygen atom, a is the distance betweeii two 
neighboring Mn atoms, k is the elastic constant.EiP 
and (73i are quantum variables corresponding to the Q2 
and Qj, modes at a manganese site i, respectively, in the 
length unit a. Thus, the ZPE per an Mn site is found 
to be i?zPE = {ft/a)\/k/mo. Since the elastic constant k 
is roughly estimated from the frequency of a phonon of 
oxygen-bond stretching as the order of 10 - 100 eV, i?zPE 
is about 0.01 - 0.04 eV, where, mo ~ 2.7 x 10-26[kg] and 
a~4A in Lai_a;Sra;Mn03. This ZPE should be com- 
pared with the energy gain from the JT coupling within 
the classical treatment of the O atoms. It is known from 
the quantum Monte Carlo calculation that the JT energy 
gain from the staggered and long-ranged distortion is the 
largest at half fUUW ^^'^ becomes smaller with increas- 
ing hole-doping £3^ E£P The JT energy gain and the ZPE 
become competing at about 10%-doping, above which 
such long-ranged distortion melts. This is why it is nec- 
essary to realize short-ranged correlations of distortions 
in the calculations. Note here that the motion of the lat- 
tice has a time scale slower than that of the electrons. If 
one is-jOn the viewpoint of the anti-adiabatic approxima- 
tion^^ interactions between the electrons and the lattice 
become instantaneous so that the effect of the lattice is 
renormalized into the Coulomb interaction, resulting in 
a simple complement to the physics of the Mott tran- 
sition. On the other hand, the adiabatic treatment of 
the lattice can consider the difference between the two 
time scales as a limiting case of the retardation effect 
in the electron-phonon interaction. A treatment of the 
lattice fluctuation from this side of viewpoint would cap- 
ture at least an aspect overlooked within the treatment 
of the Coulomb interaction only. In this paper, we treat 
the lattice classically and perform calculations simulat- 
ing lattice fluctuation by taking Monte Carlo (MC) sam- 
pling at a given energy corresponding to the above ZPE. 

Now, we examine a model Ti. = Hd + TijT + ^ph, 
where the first, second and third terms are the elec- 
tronic part, the JT interaction between the lattice dis- 
tortion and the electron, and the elastic term of the 
displacement of the oxygens, respectively. Tid is de- 
rived from the DE model with the e^-orbital degener- 
acy under the strong Hund's-rule coupling and perfect 
spin polarization as J2^jJ2^„> f^f 4^Cj„' + UY^^nnna 
with the hopping integral tj'^^ and the effective interor- 
bital Coulomb repulsion U . The orbitals d^2_y2 and 
dsz^-r^ correspond to = 1 and 2, respectively. We con- 
sider the nearest- neighbor (NN) hopping tjj = — 3/4to, 
if| -l/4io, tjf = t^l = -(+)V3/4to along the x{y)- 
direction. The hopping amplitude to is modified by the 
displacement Uij of the O atom between the NN Mn 
atoms from the center according to the Harrison's law, 
that is to = to{l — 4:ufj/a'^)'^'^/^. The most important 
effect of three dimensionality is that the hopping in z di- 



rection favors OP into the ^322^^2 orbital. This OP oc- 
curs because the transfer amplitude in the z direction is 
larger for the d^z^-r^ orbital than dj.2_y2. To mimic this 
effect in 2D lattice, we add a lower chemical potential to 
d-^^2_j.2 orbital thaiL (i2.2_j^2 so that a 3D isotropy of the 
OP is reproduced .E^P The added chemical potential can 
be viewed as the self-hopping term with an amplitude 
to for a single-layer system under the periodic boundary 
condition in z direction. The JT interaction TYjt is given 
by -g J2iiQ3iTl +q2iT.n where Ti is a pseudo-spin oper- 
ator for orbitals. q2i (qsi) is an amplitude corresponding 
to the Q2 (Qa) normal mode in the distortions of the 
octahedra formed by the oxygens surrounding i-th Mn 
atom. The elastic term iJph is k ufj, with j run- 

ning over only NN sites of i-th site. 

The method employed in this work is the Lanczos ex- 
act diagonalization of the finite-size cluster of the above 
Hamiltonian. The system size is a/To x VlO- The en- 
ergy of the electronic part is minimized by choosing 
the boundary condition, where phase shift between two 
boundaries is optimized by introducing a flux. Note that 
such boundary condition successfully reduces the finite- 
size effect .ES In fact, one can see only small differences 
betweeri results of \/T0 x \/T0 and 4 x 4 in the previous 
studyta) without JT distortions. To realize the fluctuat- 
ing lattice distortions, we combine the above exact diago- 
nalization with the MC sampling. At each update, every 
oxygen is moved from its latest position to a place which 
is uniformly random in a range with a certain width. 
It provides a possible new distortion pattern. Here, the 
width is adjusted as the acceptance ratio becomes about 
70%. Next, an exact diagonalization is performed for 
the new distortion pattern. The obtained total energy 
E'q is used to judge by the Metropolis algorithm of com- 
parison with the Boltzmann weight cxp^—Eq/ Eo) deter- 
mined from a given energy Eo whether the new pattern 
is accepted or not. We interpret that Eq characterizes 
the quantum zero-point energy. The updated pattern 
is used in the next MC step of the oxygen displace- 
ments. At each diagonalization step, (7{lu) is obtained 
from [ax{uj)+ay{u!)]/2, where ^^(cli) = 2'Ke'^Da5{uj) -\- 

iri:^(^o) 'gj^:'lf '^(^ - Kn + E'^), whose incoherent 
part is calculated by the continued-fraction-expansion 
method.Ea' Here, Da is the Drude weight, ja is a cur- 
rent operator along a-direction (a=x, y) given by ja = 

E ^M+5„ i.4v^i+Sc.y'~c\+s„y(^^^)^ ^nd |m) denotes an 
eigenstate of the system with the energy eigenvalue 
|m=0) represents the ground state. Note that the kinetic 
energy per site is given by ^ a{uj)duj. After repeat- 
ing the 200 warm-up cycles starting from the mean-field 
solution, we calculate the MC average of 17(1^) by 1800 
sampHngs. We have performed calculations at g/to=10 
and fc/io=100. These parameters are known to repro- 
duce welLtha. experimental results at half filling around 
U/to -5.bN) We here calculate ^(ajlfor U/to=6 and 
the hole-doping concentration S is 0.2,t3' which are close 
values to the ones suggested from the experimental and 
previous theoretical works. In the above parameters, to 
is found to be ~0.2 when the mean-field solution is at 
fe,g3^) ~ (0,0.031). 
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In Fig. 1, the result of the distributions in the ((72, qs) 
plane is shown. One can see that when the given energy 
in the MC procedure is small {Eo/to=0.02), fluctuations 
of distortions are restricted only near the mean-field solu- 
tion. On the other hand, when the given energy becomes 
larger {Eo/to—0.1), the distortions distribute in a larger 
area. This area includes (0,0) corresponding to the case 
where there are neither Q2- nor Qs-modc distortions. 
The system with the mean-field distortion has an equiv- 
alent point of a local energy minimum at a negative q^i 
with (?2i=0. In our simulation for i?o/io=0.1, distortions 
in the area between these points are densely produced, 
which results in the realization of a large fluctuation of 
the lattice. 

Incoherent parts of (j{uj) for two cases of the large 
and small fluctuations of the lattice together with the 
mean-field case are shown in Fig. 2. In each case of 
Eo/to — 0.1, 0.02 and (mean-field case), the weight 
has a minimum at about LUc/to ^ 5, above which the 
structure of the response to the upper Hubbard band 
due to the strong correlation appears. Incoherence be- 
low LUc corresponds to the lowest-temperature structure 
below about leV discussed in Ref. 8. Thus, to is esti- 
mated as 0.2 eV, which has an order consistent to the 
experimental indications.E9 Note that with this estima- 
tion of to, Eo for the case of Fig. 2(b) is ^--^0.02 eV, which 
is in the range of the ZPE in the above argument. For 
Eo/to = 0.02 in Fig. 2(a), the structure of a{uj) is simi- 
lar to that in the mean-field case with only small differ- 
ences in the very low-energy region. The Drude weight 
is D/to 0.036 ± 0.001 and the ratio D/N^s ^ 0.35, 
where N^g is the effective carrier density defined as 
^ Io° c(<^)dw. The differences of these values from 
the mean-field case, D/to ^ 0.04 and D/N^ff ~ 0.4 are 
small. This indicates that small fluctuations do not fully 
reproduce the situation of the strong incoherence domi- 
nant in the charge dynamics as observed in the experi- 
ments. On the other hand, in the case of Eo/to = 0.1 in 
Fig. 2(b), the Drude weight is D/to ^ 0.019 ±0.004 and 
the ratio D/Ncs ~ 0.13. This ratio is consistent with 
the experimental indication that the incoherence domi- 
nates the charge dynamics. Let us compare this Drude 
weight with the Drude weight for the free-electron sys- 
tem -DpE = ?i^n/(2TOa^), where n is the electron density 
per unit cell. The ratio Dp^/D is estimated as 40 - 60, 
which agrees with 50 - 80 of the mass enhancement in 
the experiments, w-dependence of the incoherence below 
ujc shows a leaned shape on the low-energy region and 
thus agrees with experimental indications. Increasing 
temperature changes this shape of the incoherent struc- 
ture in experiments, which is a future problem because 
the present method is based on the Lanczos algorithm. 
Note in addition that the gap-like structure in the en- 
ergy region lower than to /to '--^ 1 is seen in the mean-field 
case. A similar structure is considered to originate from 
the orbital-excitaticoLgap due to the OP induced by the 
two dimensionality EZP Our calculations for the purely 
2D case with to ~ show that the OP and the gap- 
like structure survive even when the lattice fluctuation 
is given from Eo/to = O-l- In the case of Fig. 2(b), on 
the other hand, the large lattice fluctuation forces the 



depolarized orbital, which results in the disappearance 
of the structure. 

This strong incoherence occurs as a result of the co- 
operative effect of the orbital degeneracy, the electron 
correlation and the largely fluctuating JT distortions. 
To see the importance of the complex fluctuations, let 
us compare this incoherence with cases of various ampli- 
tudes of the uniform Qa distortion including the mean- 
fleld solution (see Fig. 3). Among the cases of uniform 
Q3 distortion, the charge dynamics is the most incoher- 
ent around the energy-optimization points for the mean- 
fleld solutions. In the larger amplitude of the uniform 
Qs distortion, the orbitals become polarized and either 
the orbital d^^2_j.2 or d^2_y2 is favored. In the orbitally 
polarized state, electrons tend to move only within the 
favored orbital. This is why the charge dynamics re- 
covers the coherence for the larger amplitude of uniform 
Q3 distortion than the mean-field case. Therefore, it is 
understood that the strong incoherence obtained in the 
MC calculation for Eo/to = 0.1 cannot be reached by 
only the uniform distortion with any amplitude and 
that this strong incoherence is not realized only from a 
simple mixture of the cases with various amplitudes of 
uniform Qj, distortions, one of which gives the mean-field 
solution with the largest JT energy gain. 

Finally, two controversial interpretations about the 
small Drude weight .in experiments should be mentioned. 
Takenaka et alr^c3 made the measurement on cleaved 
surfaces and claimed that the coherence is not so small 
from the viewpoint of the obtained structure leaned on 
the low-energy region. From the measurement by Okuda 
et al& after Takenakas' work, a{uj) on the surface with 
annealing procedures removing the residual stress seems 
very similar to Takenakas' result. Note that all the 
experimental results show the common feature of this 
leaned shape up to '--^leV, which is reproduced here in 
Fig. 2(b). However, Takenaka et al. did not estimate 
the true Drude response sharply centered at uj=Q and 
instead assigned the broad incoherent response as the 
Drude weight, while the estimation of the Drude weight 
by Okuda et al. supports its smallness. Although Take- 
naka et al. attribute the controversy to the surface effect, 
our calculation suggests that pure two dimensionality re- 
covers the coherence due to the increasing OP. 

In summary, we have studied the charge dynamics in 
the perovskite Mn oxides within the framework of the 
generalized double-exchange model including the three 
elements, namely, the orbital degeneracy, the electron 
correlation, and the Jahn- Teller coupling with largely 
fluctuating distortions. As the cooperative effects of all 
the elements, the optical conductivity, whose incoherence 
dominates the charge dynamics, has been obtained. The 
result is favorably compared with the incoherent charge 
dynamics observed in the experiments. Thereby it is 
found that not only the strong electron correlation in the 
orbital system but also the Jahn- Teller interactions with 
large fluctuations are indispensable in the minimal model 
describing the properties of the compounds. The man- 
ganite provides a suitable example where the interplay 
of all these elements is crucial not only in the undoped 
Mott insulator and the phase separated system in the 
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very slightly doped caseSlli) but also in the doped metal 29) K. Takcnaka, K. lida, Y. Sawaki. S. Sugai. Y. Moritomo and 

near the Mott insulator. Contrary to the simple problem A. Nakamura: preprint |:ond-mat/99053lo| . 

of random potential without electron correlations, elec- 
tron motion has already been strongly obstructed by the 
correlation, under which the incoherence is triggered by 
the lattice fluctuation. Our results here give insights on 
such coupled effects under entanglement. In this work, 
the lattice fluctuation is treated classically and simulated 
by making such situations in the Monte Carlo sampling. 
As a future problem, direct treatment of the quantum 
nature of the lattice is desired. 
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